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1 Abstract 

The purpose of this paper is to discuss the use of 
Computer-Aided Design (CAD) geometry in a Multi- 
Disciplinary Design Optimization (MDO) environ- 
ment. Two techniques are presented to facilitate the 
use of CAD geometry by different disciplines, such as 
Computational Fluid Dynamics (CFD) and Compu- 
tational Structural Mechanics (CSM). One method is 
to transfer the load from a CFD grid to a CSM grid. 
The second method is to update the CAD geometry 
for CSM deflection. 


2 Introduction 

The process of aircraft design can be broken into 
three phasesfl]: (1) conceptual design, (2) prelimi- 
nary design, and (3) detail design. The conceptual 
design process focuses on the basic design optimiza- 
tion of features, such as weights, sizes, and overall 
performance. During the preliminary design, the fo- 
cus is on the mathematical modeling of the outside 
skin of an aircraft with sufficient accuracy. After this 
phase, the geometry is frozen, and any change could 
be costly. Detail design concentrates on the actual 
design of pieces to be fabricated. 

Often an aircraft is represented by a simple model 
during the conceptual and preliminary designs. Be- 
cause simple models are neither accurate nor com- 
plete, optimization of these models could lead to an 
impractical design [2, 3]. This shortcoming can be 
alleviated by using a high fidelity model, and the in- 
teraction among various disciplines must be modeled 
accurately. These interactions are very complicated 
and important piece of MDO. 

The strong interactions of CSM and CFD are very 
common in an MDO environment. Such interactions 
can prompt physically important phenomena such as 
those occurring in aircraft due to aeroelasticity. Cor- 
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reef modeling of these complex aeroelastic phenom- 
ena requires direct coupling of CSM and CFD. Dur- 
ing optimization of a flexible structure (e.g., wing), 
the geometry changes due to the aeroelastic effect. 
All disciplines share the same geometry, and must be 
able to consistently communicate and share informa- 
tion (e.g. deflection and load). The geometry repre- 
sentation for MDO must be accurate and suitable for 
grid generation for various disciplines such as CFD 
and CSM. To further complicate the process, model- 
ing of complicated geometrical models requires use of 
CAD systems. 

The interactions among various disciplines require 
the manipulation of the original CAD geometry that 
is stored as a set of NonUniform Rational B-Splines 
(NURBS). This paper describes two techniques to 
manipulate the NURBS geometry. In the next fol- 
lowing sections there are brief discussions on NURBS, 
parameterization of aircraft geometry, NURBS-based 
optimization, load transfer, deflection transfer, re- 
sults, and conclusions. 


3 NURBS 

This section contains a brief overview of the NURBS, 
and readers should consult [4] for a detailed discus- 
sion. A NURBS curve, R(u), can be represented as 

- , = e : =0 Nj^WiPj 
f e : = 0 N hP ( U )w t • 

The parameter, u, is bounded by < u < u max . 

The Pi are the control points (forming a control poly- 
gon), and Wi are the weights. The A/ p are the p-th 
degree B-spline basis functions defined on the non- 
periodic and nonuniform knot vector ( u ) 

U £ [u m in — Uo C Ui ^ . . .Uj C . . . ^ Uk — U max ] : 

( 2 ) 

where k is the number of knots. This completes 
the mapping between the one-dimensional parameter 
space, u, and the three-dimensional Euclidean space, 
R. A NURBS curve has five important properties: 

• It is invariant under linear transformation. 
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• A NURBS curve of order p, having no multiple 
interior knots, is p - 2 differential able. 

• The approximation is local in nature. 

• A NURBS curve is contained in the convex hull 
of its control points. 

• The NURBS approximation is variation dimin- 
ishing. 


To evaluate, the three-dimensional curve NURBS is 
commonly represented in homogeneous form as 

R(X,Y,Z) ^ R w {WX,WY,WZ,W). (3) 

So, the NURBS curve can conveniently be defined as 
a perspective map of its nonrational counterpart in 
four-dimensional space as 

n 

R w = J2 N ^ P ^ W ’ ( 4 ) 

2 = 0 

where is defined as 

P™ = {XiWi^YiWi^ZiWi^Wi}. (5) 

The basis functions can be efficiently computed by 
using DeBoor algorithm [5] as 




1 if Ui < u < Ui + i,Ui < u i+ 1 , 
0 otherwise 


( 6 ) 


Ni, p (u) = Li(u)Ni p -i(u) + Mi{u)N i+ i tP -i{u), (7) 


where 


Li{u) 


U — Ui 

u i-\-p 


= u * +p+1 ~ u 
It is agreed that 0/0 = 0. 


( 8 ) 

(9) 


A NURBS surface is a parametric surface and is 
defined as a function of two parameters as 


R(U) = {X(U),Y(U),Z(U)} t , 

U = {u 1 ,u 2 } t E [(a,b),(c,d)l 

where the components of vector, U , are the sur- 
face parameters and have no geometrical significance. 
However, for a constant u 2 , as u\ increases the point, 
R(U) moves always from one side of the surface to 
the other side. The NURBS definition for the surface 
is defined as 


R(U) 


J2i = 0 J2j= 0 Ni t p(ui)Nj t q(u2)WijPiJ 
J2i = o J2j = o 


• ( 11 ) 


where Pi j are control points (forming a control sur- 
face), Wij are the weights, and JV^ p and Nj >q are 
the p-th and q-th degree B-spline basis functions de- 
fined on the non-periodic and nonuniform knot vec- 
tor. The evaluation process is very similar to the 
NURBS curve evaluation. 


4 Aircraft Parameterization 

To use complex shapes in an MDO environment, the 
parameterization and geometry modeling must be 
compatible with existing CAD systems, and it must 
be adaptable to CFD (i.e. , block-structured or un- 
structured grids) and CSM. The four approaches for 
parameterization of an aircraft geometry are based 
on: analytical, semi-analytical, discrete, and CAD 
representation. 

Analytical methods for optimization have been 
used for a long time. These methods converts a set 
of design variables (e.g., wing sweep, thickness ra- 
tios) into a set of surfaces. Then, these surfaces can 
be used to analyze and compute the objective func- 
tion (e.g., [6, 7, 8]). This approach is very simple and 
requires a few design variables. The geometry based 
on this approach is very smooth. On the other hand, 
the geometry can take a limited form, and it is hard 
to optimize existing and free-form geometry. Inter- 
actions among disciplines are very difficult to model. 
Blair and Reich [8] have implemented a Virtual De- 
sign Process (VDM) that is integrated with full as- 
sociativity within Pro/Engineer CAD/CAM software 


In the second approach, semi-analytical, a set of 
points can describe the initial geometry, and a poly- 
nomial can model the perturbation of geometry [10]. 
Then, the coefficients of this polynomial are used as a 
set of design variables. Again, this method is simple, 
and it allows the designer to use existing geometry. 
This approach requires a few design variables, and the 
smoothness of the geometry depends on the baseline 
geometry. This approach is very difficult to general- 
ize and use in an MDO environment for a complex 
geometry. 

The third approach, discrete, is based on a discrete 
representation of the geometry. The baseline geom- 
etry creates the grids, and the position of each grid 
point becomes a design variable for the optimizer. 
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This is very easy to implement, and the geometry 
changes don’t have a limited form. The latter could 
create a problem in which the optimum design may 
not be practical to manufacture. The number of de- 
sign variables often becomes very expensive which, 
leads to high costs and a difficult optimization prob- 
lem to solve. Also it is difficult to maintain a smooth 
geometry. 

The fourth approach, CAD, is based on the op- 
timization of a NURBS representation of geometry 
(e.g., [11]). The large number of design variables of- 
ten becomes very expensive which leads to high costs 
and a difficult optimization problem to solve. The 
geometry continuity and smoothness are guaranteed. 
Also, the geometry can change locally without affect- 
ing everything else. This type of parameterization is 
flexible enough to represent a wide range of geome- 
tries. Existing complicated CAD models can be used 
as the baseline model, but modeling the interaction 
among disciplines is very difficult. 

These four approaches are summarized and listed 
in Table 1. In the next section, a NURBS-based opti- 
mization that is a hybrid approach based on first and 
fourth approaches is discussed. 


5 NURBS-Based 
Optimization 

CAD systems have been developed very rapidly and 
integrated into the design process(e.g., see [8]). Use of 
CAD systems for geometry modeling in an MDO en- 
vironment could potentially save development time. 
However, there are two drawbacks: (1) initial invest- 
ment (software and training), and (2) inability to 
calculate analytical sensitivity. The geometry rep- 
resentation in these systems is complicated due to 
a large number of entity representations. In a tra- 
ditional CAD system, the geometry is represented 
as one of many possible mathematical forms such as 
Bezier, Coons patch, B-Spline, surface of revolution, 
etc. However, one can use NURBS equations to rep- 
resent most parametric and implicit surfaces without 
loss of accuracy [4]. NURBS can represent quadric 
primitives (e.g., cylinder, cones), as well as free form 
geometry [4]. There are some CAD surfaces (e.g., 
helix and helicoidal[12]) that cannot be directly con- 
verted to NURBS. These surfaces are not common in 
most CAD system. The NURBS representation must 
be used in such as a way that it should be possible to 
altered it automatically to accommodate the changes 
due to: (1) design variable changes, or (2) structural 


and control surface deflections. 

Calculating of the sensitivity of geometry with re- 
spect to the design variables could prove to be very 
difficult. In some instances, it is possible to relate 
the NURBS control points to the design variables. 
Then the analytical sensitivity can be calculated out- 
side of the CAD system. Another way to calculate 
the sensitivity is to use finite difference, as long as 
the perturbed geometry has the same topology as the 
unperturbed one. Both methods, the analytical and 
finite difference, have their pitfalls and limitations. 

Implementation described here is based on the 
Framework for Interdisciplinary Design Optimization 
(FIDO) [13, 14] developed at NASA Langley Re- 
search Center (LaRC). The process for geometry cre- 
ation, integration, and manipulation are designed 
around the NURBS representation of the complete 
geometry. To embed this process into an optimiza- 
tion process such as FIDO, the model (High-Speed 
Civil Transport) has to be parameterized with a set 
design variables (DVs). 

For each optimization cycle, the run starts with 
a definition for a set of design variables (see Figure 
1), to build a NURBS-based geometry. Pro-Engineer 
([9]) is proposed for this step as the geometry builder. 
Once the geometry is built, the NURBS geometry 
will be deposited into a NURBS database, which will 
be shared among different disciplines. This database 
will be maintained as the baseline geometry changes. 
During the optimization process, CFD and CSM dis- 
ciplines will need the complete geometry definition. 

For each iteration of the aeroelastic loop, the CAD 
geometry is used to create CFD and CSM grids. This 
requires that both disciplines use and modify the 
same unique geometry in the geometry database. In 
each loop, the CFD grid is used to compute the aero- 
dynamic load, which is transferred to the CSM grid. 
This load will be converted to a NURBS definition 
and deposited into the database. At this point, the 
CSM grid will be generated based on the NURBS 
database. The NURBS CFD load will be mapped to 
the CSM grid, which along with the CSM grid is used 
to compute the resulting deflection. The next criti- 
cal step is to modify the original CAD geometry to 
accommodate and reflect the aeroelastic deflection. 

Each optimization loop can be summarized as 
shown in Figure 1: 

• 1. The CAD system, Pro/Engineer, will convert 
a set of design variables (D) into a NURBS geom- 
etry (G). Because the source for Pro/Engineer 
code is not available, the sensitivity has to be 
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calculated outside of the system. 

• 2. At this point, the geometry, G, and its 
sensitivity, GD, will be stored in the NURBS 
database. 

• 3. The Coordinates and Sensitivity Calcula- 
tor for Multidisciplinary Design Optimization 
(CSCMDO) ([15]) will calculate the CFD mesh 
and its sensitivity M, MD. 

• 4. The CFD code, with the capability to calcu- 
late the sensitivity, will obtain the load, which 
will be stored in the NURBS database. 

• 5. Similarly, the CSM grid generator can cre- 
ate the finite element grid and its sensitivity. 
The load stored in the NURBS database will be 
mapped onto the finite element model. 

• 6. The CSM code, with the capability to cal- 
culate the sensitivity, will obtain the deflection 
which will be stored in the NURBS database. 

• 2’. The final step is to modify the NURBS 
database to accommodate the new structural de- 
flection. 

The two steps that are defined and discussed in 
this paper in detail are: (1) converting CFD load 
to a NURBS definition, and (2) modify the NURBS 
geometry to accommodate the deflection. 


distance, d, between the r and R(w) is minimal and 
w is constrained to £ [(a,b), (c, d)]. The distance, d, 
can be written in terms of parameters w as 

d 2 (w) = f{w) = \R(w) — r\ • | R(w) — r\. (12) 

The next step is to fitting a single-value three- 
dimensional surface, F = F(U, V). This surface 
can be fitted based on a least-squares approximation 
[17, 18] that minimizes the approximation error. A 
three-dimensional curve is used as an example of the 
least-squares fitting. 

A set of points in three space, r(u), can be fitted 
by a B-Spline curve, R(u). The B-Spline equation 
can be expressed at each parameter, u, as 

?j = Ni t k( u j)Ri + N2 t k( u j ) 7?2 + ••• + N nt k(uj )R n 

(13) 

The above equation can be expressed as 

[7] = [N][R] (14) 

where 


[?T 5 •••? ^ jmax\ 

(15) 

- [i?i, i?2 , •••) Rn] T 

(16) 


6 Load Transfer 

To transfer the load to CSM elements, one needs to 
be able to compute the load at any point on the sur- 
face. Initially, the load, F(R)cfd , is computed on 
the CFD grid, ( Rcfd = {A, Y, Z} T ), which could be 
a set of structured quadrangles or unstructured tri- 
angles with the appropriate connectivity. There are 
two basic problems in fitting this data with NURBS: 

(1) the data has four dimensions, (X,Y, Z, F), and 

(2) the CFD grid could be an unstructured grid. 

The first problem can be solved by mapping the 
CFD grid to the original NURBS surfaces, hence, 
reducing the dimension from four, (X,Y, Z, F), to 
three, ( U,V,F ). The U,V are the parametric coor- 
dinates of the original NURBS surface. This infor- 
mation may be available from a CFD grid generation 
process. If not, the CFD grid points can be projected 
onto the original NURBS surface [16]. The process 
of projecting a point, r = { X , Y, Z} T , on a surface, 
R(U), can be performed by finding a w such that the 


[N] 


Ni ik (ui) 
Ni t k(uj) 
Ni,k( 'Ujmax) 


X n> k( u l) 
X n ,k(, u j) 

X n } k (^U max 


(17) 


where jmax is the maximum number of data points, 
?j. R n is the nth B-Spline control points. If jmax = 
n, the matrix [Y] is a square matrix and the control 
points can be calculated directly by matrix inversion, 

[i?] = [Y]- 1 [7]. (18) 


In this case, the resulting B-Spline curve passes 
through each data point. However, if the number 
of data points, jmax, is greater than the number of 
control points, n, the problem is over-specified. A 
least-square method can solve the problem as, 

[i?] = [[Y] T [Y]]- 1 [Y] T [7]. (19) 


The least-squares approximation for surfaces is 
very similar to the least-squares approximation for 
the curves. The minimization error can be written as 
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8 Results and Discussions 


Error = J2o E n, 

N, tP (U)N Jtq (V)W, tJ F , tJ ~\ 2 
N,, r (U)N ] , q (V)W,, ] 

( 20 ) 

where Fi j are control points for the NURBS sur- 
face representing the CFD load, Wij are the weights, 
and Ni p and Nj q are the p-th and q-th degree B- 
spline basis functions dehned on the non-periodic and 
nonuniform knot vector. N is the number of points 
in the CSM grid. This forms a system of linear 
equations that can be solved for control points of a 
NURBS surface representing the load. 


E n = 


F(U n ) - 


£' £ 


e: e; 


7 Deflection Transfer 

As described in the previous section, the load is de- 
fined on the CFD grid, Rcfd = {X,Y, Z} T . On the 
other hand, the deflection, ARcsm, i s dehned on the 
CSM grid, Rcsm = {X, Y, Z} T , which is represented 
by a set of polygons (e.g., triangles and quadrangles) 
with appropriate connectivity. The goal is to modify 
the CAD geometry definition, R(U), such that it re- 
flects the deflection produced by CSM. The algorithm 
for deflection transfer has four steps: 

1. Project the Rcsm onto the original NURBS sur- 
face. 

2. Create a NURBS surface based on the deflection, 
A Rcsm(U, V), which has the same degree as the 
original NURBS surface. 

3. Add/remove knots from the new surface to make 
it compatible with the original NURBS surface. 

4. Add the control points to the original NURBS 
surface to form the new surface. 

This algorithm has following properties: 

1. As A Rcsm approaches zero, the method will re- 
produce the original NURBS surface, R(U). 

2. Smoothness is controlled on the resulting 
NURBS surface. 

3. The results surface is a NURBS surface with the 
same degree as the original NURBS surface. 

4. It is possible to maintain the same knot vector 
as the original NURBS surface. 


The results are presented for the load and deflection 
transfers. For the load transfer, a generic pathfinder 
geometry (see Figure 2) is used. The geometry is 
a single cubic NURBS surface with 53 by 24 control 
points. There are two test cases for the load transfer. 
For the first test case, a Sin function is dehned over 
a triangular mesh (see Figure 3), which covers the 
surface of the pathfinder. The data is fitted with a 
NURBS surface. The original and interpolated con- 
tours are shown in Figure 4. The Root Mean Squares 
(RMS) error for interpolation is less than two per- 
cent, and the resulting cubic NURBS surface has 15 
by 15 control points. For the second test case, the 
pressure distribution on the surface is fitted with a 
cubic NURBS surface (see Figure 5). The resulting 
NURBS surface has 35 by 24 control points, and the 
RMS error is one percent. Figure 6 shows a close-up 
view near the inboard leading edge. 

The next two test cases demonstrate the results 
of deflection transfer. A generic High-Speed Civil 
Transport (HSCT) geometry is used for both test 
cases. This geometry is made of three surfaces: 
fuselage, inboard wing, and outboard wing. Figure 
7 shows the original NURBS surface, the deflected 
CSM grid, and the deflected NURBS surface. To 
test the limits of the deflection-transfer algorithm, 
the CSM grid has a large and unrealistic deflection. 
Figure 8 shows a side view of the same test case. 
For the last test case, the HSCT model was deflected 
by using a Sin function. Figure 9 shows the original 
NURBS surface, the deflected CSM grid, and the de- 
flected NURBS surface. Figure 10 shows a side view 
of the same test case. The resulting NURBS surfaces 
have an RMS error of less than one percent. For the 
last two test cases, even though the inboard and out- 
board were fitted separately, no gaps or overlaps were 
observed. But this could be a potential problem. 

Using the techniques described here, it is possible 
to include CAD geometry within an MDO environ- 
ment with strong interdisciplinary interactions. The 
potential drawbacks for these algorithms are: (1) cal- 
culating sensitivity, and (2) potential gaps and over- 
laps among multiple surfaces. Both need to be stud- 
ied further. 


9 Acknowledgments 

The author would like to thank Jill Vaden for proof- 
reading this paper on such a short notice. 


6 



References 


[1] Raymer, D.P. “ Aircraft Design: A Conceptual 
Approach,” AIAA, AIAA Educational Series. 
Washington D.C., 1989. 

[2] Aidala, P.V. Davis, W.H. Mason, W.H., “Smart 
Aerodynamic Optimization,” AIAA Paper 83- 
1863, 1983. 

[3] Hutchison, M.G., Huang, X., Mason, W.H., 
Haftka, R.T., and Grossman, B. “Variable- 
Complexity Aerodynamic-Structural Design of a 
High-Speed Civil Transport Wing,” AIAA-92- 
4695, September 1992. 

[4] Farin, G. “Curves and Surfaces for Computer 
Aided Geometric Design, ” Academic Press, 1990. 

[5] DeBoor, C. “On Calculating with B-Splines,” 
Journal of Approximation Theory, no. 6, pp. 50- 
72, 1972. 

[6] Bloor, M.I.G., Wilson, M.J. “Efficient Parame- 
terization of Generic Aircraft Geometry,” Journal 
of Aircraft, vol. 32, no. 6, pp. 1269-1275, Decem- 
ber 1995. 

[7] Smith, R. E., Bloor M.I.G., Wilson, M. J., 
and Thomas, A.T. “Rapid Airplane Paramet- 
ric Input Design (RAPID),” Proceedings of 12th 
AIAA Computational Fluid Dynamics Confer- 
ence, AIAA 95-1687, San Diego, CA, June 1995. 
(ftp://techreports.larc.nasa.gov/pub/techreports/ 
larc/95/NASA-aiaa-95-1687.ps.Z) 

[8] Blair, M., Reich, G. “A Demonstration of 
CAD/CAM/CAE in a Fully Associative 
Aerospace Design Environment, ” AIAA Paper 
96-1630, April 1996. 

[9] Pro/Engineer is a licensed product of Parametric 
Technology Corporation, (http://www.ptc.com) 

[10] Hicks, R.M., Henne, P.A. “Wing Design by Nu- 
merical Optimization,” Journal of Aircraft, vol. 
15, no. 7, July 1978. 

[11] Braibant, V., Fleury, C. “Shape Optimal Design 
Using B-Sphne Design, ” Computer Methods in 
Applied Mechanics and Engineering, vol. 44, pp. 
247-267, 1984. 

[12] Letcher, J.S., Shook, M. “NURBS Considered 
Harmful for Griddmg (Alternative Offered),” 4th 
International Meshing Roundtable, pp. 253-264, 
Oct. 16-17, 1995. 


[13] Townsend, J.C., Weston, R.P., and Eid- 
son, T.M. “A Programming Environment for 
Distributed Complex Computing. An Overview 
of the Framework for Interdisciplinary De- 
sign Optimization (FIDO) Project, ” NASA 
TM 109058, December 1993. (http://hpccp- 
www.larc.nasa.gov/ Hdo/homepage.html) 

[14] Sistla, R., Krishnan, R., Dovi, A. “A Computa- 
tional Framework for HSCT Design, ” the NASA 
Computational Aero-sciences Workshop, August 
13-15, Moffett Field, CA, 1996. 

[15] Jones, W. T., Samareh, J. “A Grid Generation 
System for Multi-disciplinary Design Optimiza- 
tion,” AIAA Paper 95-1689, 1995. 

[16] Samareh- Abolhassani, Jamshid, “Unstructured 
Grid on NURBS Surfaces,” AIAA-Paper 93-3454, 
1993. 

(ftp://techreports.larc.nasa.gov/pub/techreports/ 

larc/93/NASA-aiaa-93-3454.ps.Z) 

[17] Hayes, J.G., Halliday, J. “The Least-Sguares Fit- 
ting of Cubic Spline Surfaces to General Data 
Sets, ” Journal of Institute of Applied Mathemat- 
ics and Applications, no. 14, pp. 89-103, 1974. 

[18] Rogers, D.F., Fog, N. G. “Constrained B-Sphne 
Curve and Surface Fitting,” Computer-Aided De- 
sign, vol. 21, no. 10, pp. 641-648, December 1989. 


7 




Figure 1: NURBS Based Optimization 


Criteria 

Analytical 

Semi- Analytical 

Discrete 

CAD 

Design Variables 

Few 

Few 

Many 

Many 

Level of Complex Geometry 

Low 

Low 

Medium 

High 

Discipline Interaction 

Difficult 

Difficult 

Fair 

Difficult 

Smoothness of Geometry 

Easy 

Easy 

Difficult 

Easy 

Shape Limitation 

Very Limited 

Limited 

No Limits 

No Limits 

Local Control 

None 

Some 

Good 

Excellent 

Analytical Sensitivity 

Easy 

Easy 

Difficult 

Very Difficult 


Table 1: Summary of Parameterization Techniques 
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Figure 2: A Generic Pathfinder 



Figure 3: Triangular Mesh for the Generic Pathfinder 


Original Function Interpolated Function 



Figure 4: Transfer of Sin Function for a Generic Pathfinder 
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Original Function 


Interpolated Function 



Figure 5: Load Transfer for a Generic Pathfinder 


Original Function 


Interpolated Function 




Figure 6: Load Transfer for a Generic Pathfinder (Close Up) 
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Original NURBS Surface 


Original NURBS Surface 


Figure 7: Deflection of a Generic HSCT Configurations (Front View) 



Figure 8: Deflection of a Generic HSCT Configurations 
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Original NURBS Surface Original NURBS Surface 


Figure 9: Deflection of a Generic HSCT Configurations Based on a Large Deflection (Front View) 



Figure 10: Deflection of a Generic HSCT Configurations Based on a Large Deflection 
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